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Readout of the state of a superconducting qubit by homodyne detection of the output signal from a 
dispersively coupled microwave resonator is a common technique in circuit quantum electrodynamics, 
and is often claimed to be quantum non-demolition (QND) up to the same order of approximation 
as in the dispersive approximation. However, in this work we show that only in the limit of infinite 
measurement time is this protocol QND, as the formation of a dressed coherent state in the qubit- 
cavity system applies an effective rotation to the qubit state. We show how this rotation can be 
corrected by a coherent operation, leading to improved qubit initialization by measurement and 
coherent feedback. 


For most quantum information and computing proto¬ 
cols measurement is a necessary component, either to 
extract the answer to a computation, or as an operation 
in the protocol, such as for entanglement generation or 
gate operations. In addition, many protocols benefit from 
so called quantum non-demolition (QND) measurement, 
where the Hamiltonian describing the measurement oper¬ 
ator commutes with the self-Hamiltonian of the system 
[lj. As a result, perfect QND measurement maximally 
dephases the system in its eigenbasis, and the system 
state is projected onto an eigenstate when the measure¬ 
ment result is observed. Alternatively, one can think of 
a QND measurement as having only the minimal (re¬ 
quired by quantum mechanics) back action on the system 
it measures. 

In the field of circuit quantum electrodynamics 
(cQED), the state of a superconducting qubit is typi¬ 
cally measured in its eigenbasis by homodyne detection 
of the phase of the output signal through a cavity disper¬ 
sively coupled to the qubit [2, 3]. Due to the small signal 
strength exiting the cavity, it is necessary to amplify the 
signal using a low noise (near quantum limited) para¬ 
metric amplifier based on the nonlinearity induced by a 
Josephson junction [4-11]. In recent years this measure¬ 
ment scheme has been a great success, with highlights 
that include the observation of qubit quantum jumps 
[12], heralded initialization via measurement [13, 14], en¬ 
tanglement generation between qubits [15, 16], quantum 
teleportation [17], and readout fidelity greater than 99% 
[18]- 

Under the dispersive approximation, this readout 
scheme has been reported in the literature to be QND [2], 
as in the dispersive frame the qubit-cavity coupling is di¬ 
agonal and commutes with the system self-Hamiltonian. 
However, it was recently shown [19] that for a semi- 
classically driven cavity, the joint system of the qubit- 
cavity in the lab frame is an entangled state known as 
the dressed coherent state. To lowest order, this entan¬ 
glement results in a rotation of the qubit state that de¬ 
pends on the coherent state amplitude in the cavity. As 
a result, dispersive measurement as previously proposed 
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is not perfectly QND, even up to the same order of ap¬ 
proximation as the dispersive approximation, except in 
the limit of infinite measurement time. This poses prob¬ 
lems for schemes that require perfect QNDness, such as 
those performing heralded initialization or entanglement 
generation [13, 14, 20]. 

In this letter we examine the dispersive qubit read¬ 
out protocol and account for the effects of the formation 
of dressed coherent states during the protocol. In par¬ 
ticular, we describe the effective coherent qubit rotation 
that depends on both the amplitude and phase of the ap¬ 
plied cavity drive. This rotation is equivalent to a change 
of the measurement basis, and, as it is coherent, it can 
be corrected for by unitary feedback. This opens up the 
possibility for true QND measurement by introducing co¬ 
herent feedback. 

We consider a system consisting of a single qubit cou¬ 
pled to a microwave resonator (cavity), as described by 
the familiar Jaynes-Cummings Hamiltonian [21] 

H = w c a^a- ~2** z 5 ~~+ <7 + a) , (1) 

where a and a' are the usual bosonic annihilation and 
creation operators for the cavity, a z is the Pauli matrix 
whose eigenstates are the qubit logical states, are the 
qubit raising and lowering operators, w c / q are the cavity 
and qubit frequencies, g is the Jaynes-Cummings cou¬ 
pling strength, and we set H = 1 from here on. This 
Hamiltonian describes evolution in the lab frame, by 
which we mean we have not described any of the sys¬ 
tem’s evolution by a (possibly time dependent) rotation 
of Hilbert space. 

After the dispersive frame transformation, in the limit 
A = g/A < 1, where A = w q — ui c , the Jaynes-Cummings 
Hamiltonian reduces to the dispersive Hamiltonian 

Hn = io c a'a --—oy - yoya'a, (2) 

where % = g 2 / A, and we have kept terms only up to 
second order in A. The dispersive frame transformation, 
followed by the discarding of terms beyond second order 
in A is commonly called the dispersive approximation. 
Under the dispersive approximation, the system eigen- 
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states in the lab frame of equation (1) are [22] 

| g, n) = cos (Xy/n) \g, n) — sin |e, n — 1), (3) 

|e, n — 1) = cos |e,n- 1) + sin(A v / n) \g,n ), (4) 

which are referred to as the dressed eigenstates, and it 
is worth pointing out that \g, 0} = \g, 0), i.e. the dressed 
and undressed ground states are the same as \g, 0) is dark. 

The last required ingredient for dispersive qubit-state 
readout is a classical cavity drive, described by the 
Hamiltonian 


H d (t) = 2cos(wdt) (ea + e*a^) , (5) 


in the lab frame. Under the dispersive approximation 
this Hamiltonian is unaffected to lowest order in A, and 
the leading order correction term is both damped by 
the small parameter A and oscillates quickly provided 
Wd w q . Now, if we consider photons from the applied 
drive that interact with the qubit-cavity system, when 
they exit the cavity they will carry qubit information 
with them which can be used to read out the state of 
the qubit. In particular, by setting Wd = w c the qubit 
state information is contained only in the phase of the 
signal exiting the cavity, as described in Ref. [2]. These 
statements will be made more concrete shortly. 

As was shown in [19], if the system starts in either 
initial state \g/e,Q) (see appendix C for the bare excited 
state as the initial state), then after applying a cavity 
drive of the form of equation (5) for a time td the state 
of the qubit-cavity system in the lab frame will be the 
dressed coherent state g/e,a g / e (td)), defined by 


- v K/e^d)! 2 

g/e,a g / e (t d )) = e 2 


E 


Ug/e(tc l)" 

\fn\ 


\g/e,n), (6) 


where a g / e (t d) are given by 


a e (t d ) = ^ (e~ ixtd - 1) e -*("<=-x)td 

aJtd) = — (e ixtd - 1) e -H“c+x)t d 
X 

= -^sin(| (7) 


for Wd = oj c as used for dispersive readout. The phase fac¬ 
tors e -*(A c ±x)*d are due the cavity self-Hamiltonian as 
well as the dispersive interaction. The dressed coherent 
state is entangled, and correctly accounts for the corre¬ 
lations created between the qubit and the cavity during 
the applied classical drive. 

To first order in A, both dressed coherent states can be 
approximated by (see appendix A for further details) 


|p, Oig{td)) = (Iff) - \a g (td) |e)) \a g (t d )) /VM + C>(A 2 ), 

( 8 ) 


where AT = 1 + A 2 |a s / e (td)| 2 , and even for A <C 1 we can 
keep the term proportional to A|a g / e (fd)| as \a g / e {td)\ 
can be large. This approximation gives a good intuitive 
picture of the effect of an applied cavity drive on a qubit- 
cavity system described in the lab frame. The cavity is 
driven to a coherent state (as expected), while the qubit 
state is rotated a small amount. This rotation depends on 
the coherent state amplitude a g / e (td), and therefore on 
the amplitude, phase, and duration of the applied cavity 
drive. 

To connect to dispersive readout [2], we introduce a 
cavity decay mechanism via the cavity-environment cou¬ 
pling operator a + a ', described for a bare cavity by the 
quality factor Qp. For an approximately Ohmic environ¬ 
ment around the cavity frequency, such as for an open 
transmission line, the decay rate is defined in terms of 
the quality factor by k(uj c ) = w c /<5f- For a coupled 
qubit-cavity system, following the dressed decoherence 
model of [23, 24], in addition to cavity decay there will 
also be cavity-mediated qubit decay (indirect Purcell de¬ 
cay [23]). To lowest order in A (as shown in [23, 24]), for 
an approximately Ohmic environment around the qubit 
frequency, this occurs at a rate yp = A 2 (w q + x)/Q f re¬ 
gardless of the cavity photon number (photon number 
effects become relevant at higher orders of A). In exper¬ 
iment, Purcell decay can be almost completely removed 
by appropriate filtering of the cavity output at the qubit 
frequency, a technique known as Purcell filtering [18, 25]. 

In the eigenbasis of the Jaynes-Cummings Hamilto¬ 
nian, the two decay mechanisms described above amount 
to the following. “Cavity decay” is the dressed eigenstate 
transition |g/e,n-|-l) —> \g/e,n), which effectively pre¬ 
serves the qubit state, while Pu rcell d ecay is the dressed 
eigenstate transition |e,n) —> \g,n), which effectively 
preserves the photon number in the cavity. All other 
transitions have zero matrix elements with the cavity- 
environment coupling operator and are therefore forbid¬ 
den. 

Initially, let us assume that we can neglect Purcell de¬ 
cay, as would be the case if a suitable Purcell filter is 
connected to the cavity, as has been achieved in state 
of the art dispersive readout [18, 25]. In addition, we 
assume that the cavity only begins decaying after the 
state \g/e, a g / e (t d )) has been created, and we assume a 
temperature of zero. For this simplified case, given that 
only transitions of the form |c//e,+ 1) —> \g/e,n) are 
allowed, we see that the dressed coherent state will de¬ 
cay similarly to a coherent state in a bare cavity, such 
that after a time r of decay the qubit-cavity state will be 


|e, a e (td)) = (|e) + Aa*(t d ) | g)) IM^d)) /v^A7 + 0( A 2 ), 

(9) 
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cj,a g {t d )e 2 T e -' i ^-^ T ) 


- Xa g (t d )e 2 T e 


-f 


« 9 (td)e-5 T e- i ("«-^ T )/v^) + G(A 2 ), (10) 


|e,a e (£d)e 2 r e ®(“<=+x)t^ 

= f|e)+Aa e *(t d )e-f T e^+^|5) 


=(*d) 


,-f t -i{w c +x)T 


)/vm+o(\% (ii) 


which is just a dressed coherent state with a damped 
amplitude \a g / e {td)\e~ 2 T . 

Following equations (10) and (11), after a cavity decay 
of time r the initial qubit states | g/e) have been mapped 
approximately to 


While this erroneous qubit rotation scales as A and 
e~^ T and is therefore small, it is an effect that will prop¬ 
agate throughout a computation, affecting the fidelity 
of all subsequent gates and measurements. Therefore, 
the observation that it can be easily corrected for is a 
useful one, especially in architectures where qubit mea¬ 
surements are used as a means of initialization at the 
beginning of a computation [13, 14], or used during the 
computation to stabilize error correction codes [26], and 
in light of the fact that dispersive qubit readout fidelity 
approaches ever higher values [18]. 

In order to relax our previous assumption of an ini¬ 
tially closed cavity, and to consider effects beyond first 
order in A, we numerically simulate the cavity drive and 
decay process with the Purcell filtered master equation 
(see appendix B for further details, and appendix D for 
simulations with Purcell decay) 


I g) 

(| fl ) - A|a|e- i( ^ + " /2) e-^ td e-5 T e- iw ^ | e >) /V^r), 

( 12 ) 


|e) -> 

(|e> + \\a\e i ^ + 7 l/ 2 ) e iu ’* td e-% T e i < T |<?)) /yW), 

(13) 


p(t) = -i \H T {t),p(t) 

+ /c((l + n t h(wc,T))D 


ac 


+ n t h(w c , T)V 


(16) 

Pit), 


where the operator ac of equation (B9) describes only 
cavity decay, ?r t h(w,T) is the Bose distribution at fre¬ 
quency u> and temperature T, and T>(x)p is the dissipator 
defined by 


where we have used the fact that a g / e (id) differ in phase 
only to define a = |a fl / e (£d)| and the phase factor — (ip e + 
7t/2) of —ie*, as well as the shifted cavity frequencies 
u> g / e = <jJ c -F 2 an d w g/e = w c T X- The measured phases 
of the output signal jump sharply when the drive pulse 
is turned off, with the coherent states rotating around 
phase space at frequencies u g / e for times t < i d and a t 
frequencies u)' g , e for times r = £ — i d > 0 . 

The qubit-state maps of equations (12) and (13) can 
alternatively be understood as qubit-cavity interactions 
during the cavity drive changing the basis of qubit mea¬ 
surement, with measurement of the cavity frequency 
uig/uj' g for a total time i d + T corresponding to the qubit 

state |g() —A|ci!|e _ 5 T e _ ^ ¥, ' + 7 r / 2 +w s td+u ’s T ) |e) andasimilar 
result for the measurement of w e /w'. Only in the r — > oo 
limit does the basis of measurement become {| < 7 ), | e)} 
and the measurement QND. 

The key observation in this work is that the change 
of basis of measurement is a coherent rotation error ap¬ 
plied to the qubit output state. As this error in the qubit 
state is coherent, it can be actively corrected for by a sin¬ 
gle qubit rotation that applies the inverse of the unitary 
map of equation (12) or (13). This conditional rotation 
to correct the output state is defined by the unitary op¬ 
erators 

U g (£ d , r, |a|) = exp {i (cos(E g )a y - sin(E 9 )<r x ) 6} , (14) 

ll e (i d , r, |a|) = exp {i (cos(E e )<r y - sin(E e )<r x ) 8} , (15) 

where E 9 = oj g td + a>' g T + ip e + ir/2, S e = tp e + it/2 + 
w e id + w'r, and tan(0) = X\a\e~^ T . 


V(x)p = xpx* — ^ {x^x, p} . (17) 

Here the total system Hamiltonian describes the full 
Jaynes-Cummings interaction between the qubit and the 
cavity as well as the classical cavity drive, and is given 

by 

H T (t) = u> c a)a -A1<j 2 q. g (o-~a^ + a + a) 

+ (ee^ dt a + e* e - iWdt a t ) 0 (i - i d ), (18) 

where Q(x) is the Heaviside step function. We simulate 
the evolution for the initial states |.g, 0 ) and |e, 0 ), with 
the temperature set at either T = 0 or T = 100 mK. The 
results are shown in FIG. 1. 

FIG. 1(a) shows the cavity occupation during the read¬ 
out protocol. As expected, after an initial ring-up phase, 
once the drive is turned off the cavity occupation decays. 
Decay st ops o nce the steady state is reached, which is 
| g, 0) or |e, 0) for T = 0 and thermally broadened ver¬ 
sions of these states for T = 100 mK. FIG. 1(b) shows 
1 — P(i), where P(i) is the purity of the qubit state, de¬ 
fined for a reduced qubit state p(t) by P(i) = Tr[p(f) 2 ]. 
Unit purity indicates a pure state. As can be seen, for 
T = 0 the states remain very close to a pure state at 
all times, verifying the analytic results of equations ( 10 ) 
and (11). Even for T = 100 mK the states remain > 90% 
pure for either initial state. 

As the states remain mostly pure during the protocol, 
it is possible to correct the qubit state error by the uni- 
taries of equations (14) and (15), as described previously. 
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Figure 1: (a) Cavity occupation, and (b) 1 — P(t) for T = 0 and T = 100 mK for both initial states. 1 — T v (r) and 1 — (t) 
are shown in (c) for T = 0 and (d) for T = 100 mK. A drive strength of |e|/27r = 0.04 GHz, a cavity decay rate of 1 /k = 100 
ns, and |A| =0.1 were used for these simulations. 


To quantify this correction we use the overlap between 
the desired state (|g) or |e)) and the simulated reduced 
qubit state p(t). We measure the overlap before correc¬ 
tion 


Fv{t) = Tr [|zz) (u\ p(t )\, (19) 

where the subscript v € { g , e} indicates whether we 
started in | g, 0) or |e, 0), and the overlap after correction 


J, C (r)=Tr 


W) H U„p{t)Ul 


( 20 ) 


FIG. 1(c) shows the overlap error for both the uncor¬ 
rected and the corrected state for T = 0, and as can be 
seen T^(t) > T v (t) for all time (to within numerical 
precision of the simulations). For T = 100 mk, as shown 
in FIG. 1(d), this is not the case, as within roughly 75 ns 
the qubit state loses enough coherence that the unitary 
correction actually worsens the overlap. 

For both system temperatures the greatest benefit 
from correction is seen early on in the decay time, long 
before the cavity occupation has reached steady state. 
Typically one would wait for the cavity to be unoccupied 
before further operations on the qubit are preformed, as 
cavity photons are still interacting with the qubit. How¬ 
ever, in set-ups with tunable coupling between the cavity 
and the qubit [27-29], it would be possible to turn off the 
interaction between the cavity and the qubit once enough 


measurement data has been accumulated and then cor¬ 
rect the final state of the qubit. In this way one would 
could achieve both more accurate and faster initializa¬ 
tion of the qubit state via measurement and unitary cor¬ 
rection. A similar initialization scheme involving both 
cavity and qubit control has recently been implemented 
[30].' 

The analytic expressions of equations (14) and (15) for 
the correction unitaries correctly calculate the amplitude 
of the rotation, described by the angle 6. Unfortunately, 
due to higher order nonlinear effects in the full Hamil¬ 
tonian the analytic phases T, g and E e do not give good 
results. To solve this problem, we performed a brute force 
optimization over the phase of the rotation to obtain the 
excellent results shown in FIGs. 1(c) and 1(d). 

In conclusion, we have shown that during the most 
commonly used dispersive readout protocol for supercon¬ 
ducting qubits a coherent rotation error is applied to the 
qubit, and the measurement scheme is not QND for any 
finite measurement time. This coherent rotation causes 
errors in repeated measurements and in qubit initializa¬ 
tion; however, as we have shown, it can be corrected for 
by unitary feedback. This correction is most advanta¬ 
geous early on in the decay time, and in experiments with 
tunable qubit-cavity coupling our scheme shows promis¬ 
ing results for faster and more accurate qubit initializa¬ 
tion. 
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Appendix A: First Order Approximation of the 
Dressed Coherent States 

For the ground qubit dressed coherent state, we have 



-AVn|e,n—1) j +C>(A 3 ). (Al) 

Tracing out the cavity we obtain 
Tr c \g,a)(g,a\ = (l - A 2 |a| 2 ) \g) (g\ + A 2 |a| 2 |e) (e| 

- Aa* |ff) (e| - Aa |e) (ff| + £>(A 3 ) 

= (Iff) - Aa |e» «ff| - Aa* (e|) - A 2 |a| 2 1 g) (g\ + 0{ A 3 ). 

(A2) 

We now approximate this minimally mixed state by a 
pure state to obtain 

Iff, a) = (Iff) - \at\e))/VN+ 0{\ 2 ), (A3) 

where Af = 1 + A 2 |a| 2 is the normalization. 

Now for the excited qubit dressed coherent state, we 
begin with 



+ Av^TT |.g, n + 1) j + 0(X 3 ). (A4) 

Tracing out the cavity we obtain 

Tr c |e,a)(e,a| = (l - A 2 (l + |a| 2 )) |e) (e| 

+ A 2 (l + |a| 2 ) |ff) (g\ + Aa* | g) (e| + Aa |e) (g\ + 0( A 3 ) 
= (|e) + Aa* |ff)) ((e| + Aa(ff|) 

+ A 2 I. 9 ) (ff| - A 2 (l + |a| 2 ) |e) <e| + 0(A 3 ). (A5) 

Making the same approximation as for the previous case, 
we arrive at the pure state 

^) = (|e) + Aa*|ff))/^M + 0(A 2 ), (A 6 ) 

where, as before, A f is the normalization. 
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Appendix B: Numerical Master Equation 

To derive the master equations of equations (16) and 
(Dl) we begin by considering the system-bath Hamilto¬ 
nian 

H e = H T (t) + ^2 Vkblh + *22 9k (a + a f ) (b k + b^j 

k k 

(Bl) 

where the second term in equation (Bl) is the bath self- 
Hamiltonian, and the third term is the system-bath cou¬ 
pling. To derive an effective master equation for the sys¬ 
tem, it is appropriate to work in the instantaneous eigen- 
basis of however, to simply things we will derive 

the master equation in the eigenbasis of H of equation 
(1). The difference between these two is a frame transfor¬ 
mation by a time dependent cavity displacement, which 
for the parameter regime under consideration is inconse¬ 
quential to the applicability of the master equation ob¬ 
tained. 

Following the procedure of [23] we derive an effective 
evolution equation for the system (ignoring the coherent 
evolution for the time being) 

P(t) = 

C Jk C nm(\j) ( k \ p{t) \ m ) W - \™) (n\ I j) (k\ p(t)) 
j,n 
k>j 
m>n 

g iAjkt 

+ C okC* m (\j) (k\ p(t) |to) (n\ - p(t) |rn) (n\ \j) {k\ ) 

j,n 

k>j 

m>n 

ds (b{0)V {s)} e iA ™ f 
+ CjkCnm^\k) 01 p(t) \n) (rn\ - | n) (m| | k) {j\ p(t)) 

j,n 

k>j 

m>n 

k t 

+ H C* k C nm (^\k) (j| pit) | n) {m\ - p{t) \n) (m| | k) (j \) 

j,n 

k>j 

m>n 

(B2) 

where {Ij)} is the eigenbasis of H ordered in increas¬ 
ing eigenenergy and given to first order in A by equa¬ 
tions (3) and (4), Aj k is the frequency difference be¬ 
tween the j’th and fc’th eigenstate (Bohr frequency), 
b(s) = 9 kb k e~ lrikt is the time dependent bath low¬ 
ering operator, and Cj k = (j| (a + a') \k). We have also 
assumed that the bath state is a stationary state of the 
bath self-Hamiltonian. 




ds /&t(0)6(s) 


o i{^jk A„ 


ds (b\s)b{0 ) 



i(Ajk A n 


ds ib(s)P{0 ) 
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Unlike in [23], it is not possible to make a rotating 
wave approximation, as in the parameter regime under 
consideration the eigenspectrum of H has many nearly 
degenerate transitions. Instead, it is possible to derive a 
Lindblad form master equation in a way similar to that 
done in the singular coupling limit [31] by assuming that 
the nearly degenerate transitions are actually degener¬ 
ate. We notice that the coefficients Cjk are nonzero if 
I j) = |ff/e, n) and |fc) = \g/e,n± 1 ), or if \j) = \g/e,n) 
and \k ) = \e/g,n), while all other Cjk are zero. The 
former case is what we have been calling cavity decay, 
while the latter case is Purcell decay, and for each decay 
type the energy difference A jk between adjacent states 
is approximately constant (co c in the former case and w q 
in the latter). 

Therefore, we can split the sums of equation (B2) into 
two parts, and make a secular approximation to neglect 
the fast oscillating cross terms between decay types to 
arrive at the equation 

P(t) = 

c 

E CjkC* nm D | j) (k \, |to) (n\ p(t) 

j,n 

k>j 

m>n 

X (l + n t h(w c ,T)) j{u c ) 
c 

+ E C* k C nm V \k) (j\,\n) (m\ p(f)n t h(wc,T)J(w c ) 

k>j 

m>n 

V 

+ £ C 3k C* nm v[\j) (k\,\m) <■ n\\p(t ) 

j,n 

k>j 

m>n 

X (l +n t h(w q ,T)^ J(w q ) 

V 

+ E C jk C nrnV \k) {j\ ,\tl) { m \ p(t)n th (w q , T) J(w q ), 

j,n 

k>j 

m>n 

(B3) 

where we have ignored the Lamb shifts, and the super¬ 
scripts C and V indicate summation over cavity decay 
transitions and over Purcell decay transitions respec¬ 
tively. For more compact notation, we have also defined 
the two operator “dissipator” 

v [Or,0 2 ] p(t) = d lP (t)d 2 - l -{d 2 d llP (t)}. (B4) 
Following the usual procedure of Fermi’s golden rule we 


have made the identification 


j°° ds((b(s)b\0))e~ iut + 

(b(0)b\s))e^) 

= (l + n t h(w,T)) J(w), 

(B5) 

J°° ds((b!(s)K 0 ))e“ 4 + ( 

&*(())&(«)) e"***) 

= nth(w,T)J(w), 

(B6) 


where J(to) is the spectral density of the bath, and 
n t h (w, T) is the Bose function evaluated at frequency ui 
and temperature T. 

If we choose a global Ohmic spectral density J(w) = 
uj/Qf, then using the relations 

E C i k U) 

j,k>j 

E c v i fc ) 01 = E c h i fc ) oi = (B7) 

j,k>j j,k>j 

we can write equation (B3) in Lindblad form 

pit) = (l + n t h(w c ,T)) (kV d c +7 p£> dp ) p(t) 

+ n t h(w q ,T) [kD a}, + 7 pX> a\, ) p(t) (B 8 ) 

where n = u> c /Qf is the cavity decay rate and yp = 
A 2 Wq/<3F is the Purcell decay rate. The operators dc 
and dp are defined by 

C v 

ac = E C > k \3) (k\ = a~ E C j k W ^ ’ ( B9 ) 

j,k>j j,k>j 

V C 

dp = E Cjk I j) (k\ = a - E C jk | j) (k\ , (BIO) 

j,k>j j,k>j 

and describe the cavity and Purcell decay processes re¬ 

spectively. 

Equation (B 8 ) is valid for t <C where 

^me = 1 /wmax describes the timescale over which 
e *(A 3 - fc -A nm )t _ e *w max t j g no l on g er unity, with w max 
the largest degeneracy between transitions that were as¬ 
sumed to be degenerate. For cavity decay w max oc N\ A 2 , 
while for Purcell decay w max oc N\, with N the photon 
number of the largest occupied state. For the parame¬ 
ters under consideration (y < 10 MHz, A < 1CU 1 ) we 
have t me ~ 10/N /is for cavity decay, and t me oc 100 /N 
ns for Purcell decay. For typical experimental parameters 
the simulation length is on the order of 100 s of nanosec¬ 
onds, well below the limit imposed by cavity decay for 
reasonable values of N, but unfortunately on the same 
order as that set by Purcell decay. As such, we expect 
that the simulations of equation (B 8 ) are somewhat non¬ 
physical; however, as the Purcell decay rate is quite small, 
we do expect the results to remain mostly trustworthy. 

If instead of choosing a global Ohmic spectral density, 
we choose a spectral density for which J(oj c ) = u c /Qf 



but J(w q ) « 0, such as has been achieved in contempo¬ 
rary dispersive readout schemes by a filter [18, 25], then 
the resulting Lindblad equation contains only cavity de¬ 
cay, and is given by 


pit) = n(l +n t h(w c ,T)jx> a c p{t) 

+ «nth(wc,T)X> aj, p(t). 


(BH) 


This is the master equation used in the main text to 
described dispersive readout in a Purcell filtered system. 


Appendix C: Readout of the Undressed Qubit 
Excited State 

In the main text we considered starting with the ex¬ 
cited initial state |e, 0), which is the first excited state in 
the energy eigenbasis of the qubit-cavity system. How¬ 
ever, for some quantum information protocols it is ad¬ 
vantageous to work in the qubit logical basis, where the 
logical excited state is the state |e, 0). Such a state can 
be prepared by initializing the qubit-cavity system into 
its ground state \g , 0), followed by a short non-adiabatic 
qubit pulse that flips the state of the qubit. 

If we now attempt to read out the state of the qubit, 
from Ref. [19] we see that starting in the state |e, 0), after 
an applied cavity drive of the form of equation (5) the 
state of the system is 

ll'(fd)) = cos (A) |e, a e (td)) 

- e lG(td) sin (A) Ul \g) & > d D(a' g {t d )) |1), 

(Cl) 


After sufficiently long measurement of the frequency 
w e , the state of equation (Cl) will have collapsed to the 
dressed coherent state |e, a e {t d )), and the rest of the 
readout and coherent feedback protocol can occur as 
described in the main text, without further modification. 

Appendix D: Purcell Decay 


To examine the effect of Purcell decay we simulate the 
unfiltered master equation 

= H T (t),p{t ) 

+ k ((1 + n th (w c ,T))2? d c +n t h (w c ,T)2? a£ ) pit) 

+ 7p ((! + nth(w q ,T))X> d v + n th (uj q ,T)V d\, pit), 

(Dl) 


where the operator a-p of equation (B10) describes Pur¬ 
cell decay. As before, we simulate the evolution for the 
initial states |g,0) and |e, 0), at a temperature T = 0. 
The overlap error for the uncorrected and corrected states 
are shown in FIG. 2. As can be seen, for the initial state 
|< 7 ,0) the results remain unchanged, as Purcell decay only 
minimally affects the decay of states of the form \g,a). 
Coherent feedback can still be used in this case to correct 
the qubit state, with excellent improvement in the over¬ 
lap. For |e, 0) the results are quite different from that 
of FIG. 1(c), as Purcell decay drives the system to the 
global ground state | g, 0) and, as such the qubit state 
decays. In this case the coherent correction does little 
good, as the qubit state is lost by Purcell decay. 


where ag(td) = e* (e lxtd — l) /x, D(/3) is the usual 

displacement operator for a harmonic oscillator, Z7 d = 
exp (A (<f + a — } is the dispersive frame transfor¬ 

mation operator, and e lG i td> is a relative qubit phase 
whose form is unimportant [19]. From equation (Cl) we 
see that the final state contains a component for which 
the qubit is in its ground state, and the frequency of the 
cavity signal for this component will be close to u> g , mea¬ 
surement of which indicates the qubit is in its ground 
state. This introduces the possibility of misidentifying 
the qubit state; however, the amplitude of the ground 
state component is small as it as scales with sin 2 (A) ss A 2 . 
Nevertheless, this sets a fundamental limit for the read¬ 
out fidelity of the undressed excited state via standard 
dispersive readout as presented here, and partly explains 
the unequal readout fidelities reported in [18]. 



Figure 2: 1 — J- U (r) and 1 — J- G (t) for T = 0. A drive strength 
of |e|/27 t = 0.04 GHz, a cavity decay rate of 1 /k = 100 ns, 
and |A| =0.1 were used for this simulation. 


[1] V. B. Braginsky and F. Y. Khalili, Quantum Measure- [2] A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and 
ment (Cambridge University Press, New York, NY, USA, R. J. Schoelkopf, Phys. Rev. A 69, 062320 (2004). 

1992), ISBN 0-521-41928-X. [3] G. Johansson, L. Tornberg, and C. M. Wilson, Phys. Rev. 

































9 


B 74, 100504 (2006). 

[4] I. Siddiqi, R. Vijay, F. Pierre, C. M. Wilson, M. Metcalfe, 

C. Rigetti, L. Frunzio, and M. H. Devoret, Phys. Rev. 
Lett. 93, 207002 (2004). 

[5] M. A. Castellanos-Beltran and K. W. Lehnert, Appl. 
Phys. Lett. 91, 083509 (2007). 

[6] N. Bergeal, F. Schackert, M. Metcalfe, R. Vijay, V. E. 
Manucharyan, L. Frunzio, D. E. Prober, R. J. Schoelkopf, 
S. M. Girvin, and M. H. Devoret, Nature 465, 64 (2010). 

[7] L. S. Bishop, E. Ginossar, and S. M. Girvin, Phys. Rev. 
Lett. 105, 100505 (2010). 

[8] M. D. Reed, L. DiCarlo, B. R. Johnson, L. Sun, D. I. 
Schuster, L. Frunzio, and R. J. Schoelkopf, Phys. Rev. 
Lett. 105, 173601 (2010). 

[9] J. Y. Mutus, T. C. White, R. Barends, Y. Chen, 
Z. Chen, B. Chiaro, A. Dunsworth, E. Jeffrey, J. Kelly, 

A. Megrant, et al., Applied Physics Letters 104, 263513 
(2014). 

[10] K. O’Brien, C. Macklin, I. Siddiqi, and X. Zhang, Phys. 
Rev. Lett. 113, 157001 (2014). 

[11] T. C. White, J. Y. Mutus, I.-C. Hoi, R. Barends, 

B. Campbell, Y. Chen, Z. Chen, B. Chiaro, 
A. Dunsworth, E. Jeffrey, et al., Applied Physics Letters 
106, 242601 (2015). 

[12] R. Vijay, D. H. Slichter, and I. Siddiqi, Phys. Rev. Lett. 
106, 110502 (2011). 

[13] J. E. Johnson, C. Macklin, D. H. Slichter, R. Vijay, E. B. 
Weingarten, J. Clarke, and I. Siddiqi, Phys. Rev. Lett. 
109, 050506 (2012). 

[14] D. Riste, J. G. van Leeuwen, H.-S. Ku, K. W. Lehnert, 
and L. DiCarlo, Phys. Rev. Lett. 109, 050507 (2012). 

[15] D. Riste, M. Dukalski, C. A. Watson, G. de Lange, 
M. J. Tiggelman, Y. M. Blanter, K. W. Lehnert, R. N. 
Schouten, and L. DiCarlo, Nature 502, 350 (2013). 

[16] J. M. Chow, J. M. Gambetta, E. Magesan, D. W. Abra¬ 
ham, A. W. Cross, B. R. Johnson, N. A. Masluk, C. A. 
Ryan, J. A. Smolin, S. J. Srinivasan, et al., Nat Commun 
5, 4015 (2014). 


[17] L. Steffen, Y. Salathe, M. Oppliger, P. Kurpiers, M. Baur, 

C. Lang, C. Eichler, G. Puebla-Hellmann, A. Fedorov, 
and A. Wallraff, Nature 500, 319 (2013). 

[18] E. Jeffrey, D. Sank, J. Y. Mutus, T. C. White, J. Kelly, 
R. Barends, Y. Chen, Z. Chen, B. Chiaro, A. Dunsworth, 
et al., Phys. Rev. Lett. 112, 190504 (2014). 

[19] L. C. G. Govia and F. K. Wilhelm (2015), 
arXiv: 1506.04997. 

[20] C. L. Hutchison, J. M. Gambetta, A. Blais, and F. K. 
Wilhelm, Canadian Journal of Physics 87, 225 (2009). 

[21] E. Jaynes and F. Cummings, IEEE Proc. 51, 89 (1963). 

[22] S. Haroche and J.-M. Raimond, Exploring the Quantum: 
Atoms, Cavities, and Photons (Oxford University Press, 
Oxford, 2006). 

[23] F. Beaudoin, J. M. Gambetta, and A. Blais, Phys. Rev. 
A 84, 043832 (2011). 

[24] E. A. Sete, J. M. Gambetta, and A. N. Korotkov, Phys. 
Rev. B 89, 104516 (2014). 

[25] M. Reed, B. Johnson, A. Houck, L. DiCarlo, J. Chow, 

D. Schuster, L. Frunzio, and R. Schoelkopf, Appl. Phys. 
Lett. 96, 203110 (2010). 

[26] A. G. Fowler, M. Mariantoni, J. M. Martinis, and A. N. 
Cleland, Phys. Rev. A 86, 032324 (2012). 

[27] J. Wenner, Y. Yin, Y. Chen, R. Barends, B. Chiaro, 

E. Jeffrey, J. Kelly, A. Megrant, J. Y. Mutus, C. Neill, 
et al., Phys. Rev. Lett. 112, 210501 (2014). 

[28] Y. Chen, C. Neill, P. Roushan, N. Leung, M. Fang, 
R. Barends, J. Kelly, B. Campbell, Z. Chen, B. Chiaro, 
et al., Phys. Rev. Lett. 113, 220502 (2014). 

[29] S. Zeytinoglu, M. Pechal, S. Berger, A. A. Abdumalikov, 
A. Wallraff, and S. Filipp, Phys. Rev. A 91, 043846 
(2015). 

[30] K. Geerlings, Z. Leghtas, I. M. Pop, S. Shankar, L. Frun¬ 
zio, R. J. Schoelkopf, M. Mirrahimi, and M. H. Devoret, 
Phys. Rev. Lett. 110, 120501 (2013). 

[31] H.-P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford University Press, 2006). 



